A Geometric Theory of Diblock Copolymer Phases 
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We analyze the energetics of sphere-like micellar phases in diblock copolymers in terms of well- 
studied, geometric quantities for their lattices. We argue that the A15 lattice with Pm3n sym- 
metry should be favored as the blocks become more symmetric and corroborate this through a 
self-consistent field theory. Because phases with columnar or bicontinuous topologies intervene, the 
A15 phase, though metastable, is not an equilibrium phase of symmetric diblocks. We investigate 
the phase diagram of branched diblocks and find that the A15 phase is stable. 



The ability to control the self-assembly of complex lat- 
tices by manipulating molecular architecture remains an 
essential aspect in the creation of new, functional ma- 
terials. With only a few tunable parameters, diblock 
copolymer melts exhibit a wide variety of equilibrium 
phases which can be understood via the mean-field Gaus- 
sian chain model of "AB" diblock copolymers composed 
of immiscible A and B blocks P, |2|. Indeed, in a sys- 
tem where the A and B-blocks are otherwise identical, 
there are only two thermodynamic variables, 0, the vol- 
ume fraction of A type monomers, and x-^j where x 
is the Flory-Huggins parameter characterizing the repul- 
sive interactions between the A- and B-type monomers 
and N is the degree of polymerization In this let- 
ter, we present a model which predicts that the the A15 
(shown in Fig. ^) lattice of diblocks is stable relative 
to other sphere-like phases for sufficiently large (f) or, in 
other words, sufficiently symmetric diblocks. We corrob- 
orate this prediction by recalculating the phase diagram 
for symmetric diblocks (Fig. ^) via a self-consistent field 
theory (SCFT) for diblock copolymer melts 0. 

The "classical" diblock phases are well-understood: 
near the order-disorder transition (ODT), Leibler devel- 
oped a Landau-like theory in the weak-segregation regime 
to establish the stability of a body-centered cubic (BCC) 
phase, a hexagonal phase of columns and a lamellar phase 
[5j. Moreover, Semenov's picture of spherical micelles 
interacting through a disordered copolymer background 
when <C 1 accounts for the appearance of the face- 
centered cubic (FCC) lattice near the ODT in the mean- 
field phase diagram m. The more exotic gyroid phase 
was discovered 0, 11,1^1 and was explained successfully 
by Matsen and Schick via SCFT In our study of the 
A15 lattice, we find that the hexagonal and gyroid phases 
intervene and thus there should be no stable A15 lattice 
for simple diblocks. However, sphere-like topologies are 
favored by branched diblock copolymers [lOllllI ll2l | and 
dendritic polymers |l3L H^ : with this in mind we pre- 
dict that sphere-like phases are stabilized and that the 
A15 phase is a ground state for this class of structures. 
By implementing, to our knowledge, the first full SCFT 
treatment of branched molecules (shown in Fig. ^) we 
have verified our theory. 

In the dilute regime, Semenov's picture treats each mi- 



celle as an undistorted sphere so that the outer block ex- 
tends to a spherical unit cell of radius Rs ■ This unit-cell 
approximation provides a lower bound for the free en- 
ergy of the sphere-like phases However, when the 
micelles assemble into a periodic structure, the incom- 
pressibility of the melt allows no interstitial gaps and the 
spheres must deform into the Voronoi or Wigner-Seitz 
cells of the corresponding lattice. Our analysis will be in 
the strong-segregation limit, away from the ODT where 
fluctuation effects are less important. For small 0, the 
A-block will form the center of each sphere-like micelle. 
As a gedanken experiment, we increase (j) while maintain- 
ing the topology and keeping the A-block on the inside. 
In the most extreme limit, — > 1 and the A-blocks are 
surrounded by a vanishingly thin B-block coat. 

The free energy of strongly segregated diblock config- 
urations of spherical geometry arises from the tension 
in the AB interface and the stretching of the polymers. 
Though it may appear that the curvature of the inter- 
face must be considered, in neat systems the curvature 
energy is completely accounted for by stretching through 
the incompressibility constraint 0. The energy of the 
AB interface is proportional to its surface area S (and 
y^X the strong segregation limit 0, 0|). As illus- 
trated in the insets to Fig. 12 when </> grows, the A-blocks 
in the center of the micelle fill a large volume and hence 
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FIG. 1: In (a) we show the cubic unit cell of the A15 lattice. 
The black spheres form a BCC lattice and the A15 lattice 
has, in addition, dimers on each face, in grey. The diblock 
topologies we consider are (b) linear and (c) branched with 
three generations. 
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FIG. 2: A plot of the AB interfaces in the BCC (ImSm) phase 
obtained by extracting the surface where the concentration of 
each block is equal. Both figures were extracted for X-^ ~ 40. 
In (a) — 0.222 and the interface is almost spherical. In (b) 
(j) — 0.45 and the interface takes on the shape of the Voronoi 
cell. The solids lines represent the entire BCC cell, while the 
dashed lines represent the reduced Voronoi cell that would 
arise in the </> — > 1 limit. The insets depict (a) the spherical 
interface limit (small (j)) and (b) the fiat interface limit (large 



the curvature of the AB interface is small. As a result, 
the interface is easily deformed and we expect that it 
will adopt the shape of the Voronoi cell. A dimcnsionlcss 
measure of the surface energy ^/xS is 



(1) 



where Ss is the area of the interface for a spherical mi- 
celle of the same volume. This ratio may be rewritten as 
Ax = 7(-'^)/(367r)^/'^ where "f{X) depends on the lattice 
X and relates the area to the volume V of the Voronoi 
cell: E — 'y{X)V^^^. It is conjectured that an A15 lattice 
of equal volume cells minimizes 7 16] and it is known that 
7(A15} = 5.297 < 5.315 = 7(BCC) < 7(FCC) = 5.345 
|l7l llSj l . This fact can be used to argue the stability of 
the A15 phase in fuzzy colloidal systems 0| and is cru- 
cial to the analysis here. The area per unit volume scales 
as Rg^ and so the surface tension contributes A/Rs to 
the free energy density. 

The second contribution to the free energy arises from 
the stretching of the polymers to fill the cell. Due to 
incompressibility, the number of chains in each wedge 
of solid angle dO, is proportional its volume, ^R^{Q)(Kl, 
where R{Q) is the distance to the boundary. Since the 
stretching energy for both blocks of chains in this wedge 
is proportional to R^{Vl) (isf. 
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is a dimensionless measure of the stretching relative to 
the stretching in a sphere of radius Rs with the same 
volume as the Voronoi cell H. This integral is propor- 



tional to the moment: 
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which has been well-studied |20j in the context of the 
quantizing problem: G(BCC) = 0.078543 < 0.078745 = 
G(FCC) = G(A15) 113. The stretching energy per chain 
(or per unit volume) scales as R?g and so the free energy 
contribution from stretching is proportional to The 
prefactor of this term will, in general, depend on the 
topology of the micelles and (f> [l5j but for the structures 
we consider, this prefactor does not change. 

Putting these two effects together, we find a Flory-like 
free energy density: 
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Minimizing over Rs, we find that f = fo [Xj^y^^] (x 
^2/3(^1/3 ^i^QYe fo is the energy of the spherical mi- 
celle 01 . It follows, for instance, that as cj) grows, the 
BCC lattice is lower in energy than the FCC lattice since 
7(BCC) < 7(FCC) [13. We find that the A15 lattice 
is lower in energy still: Jais = 1.0707/o < 1.0722/o = 
fscc < Ifcc — l-0772/o. It should be noted that a 
more detailed, self-consistent treatment 0| yields the 
same free energy (0J and the same numerical results given 
here for the BCC and FCC lattices. We have thus estab- 
lished a fundamental and unavoidable geometric frustra- 
tion between the stretching and interfacial energies - the 
BCC lattice minimizes G and the A15 lattice minimizes 
7. We note that since the A15 lattice has two distinct 
sites, it is possible for their corresponding cells to adopt 
different volumes to lower the energy further. Indeed, by 
adjustin g th e Voronoi cells G(A15) and 7(A15) can be 
lowered jl^ I2H . ultimately lowering /^is by a further 
0.03%. 

As (j) decreases the minority region shrinks and the 
curvature grows, leading to a larger restoring force to 
the optimal, spherical shape. In this case the connec- 
tivity of the A-blocks to the B-blocks would require a 
more detailed analysis, such as the "kinked path" ansatz 
of Milner and Olmsted 0| . Since the BCC lattice mini- 
mizes the stretching energy it is not surprising that it is 
the equilibrium phase at smaller (f). 

Because we assumed that the AB interface is perfectly 
flat, our calculation puts upper limits on the free ener- 
gies of these sphere-like phases. Nonetheless, our argu- 
ment suggests that as cj) grows, the A15 phase should 
become the more stable sphere-like phase. To critically 
test our theory, we employed the SCFT implementation 
of Matsen and Schick to compare the free energy of a 
phase with the PmSn symmetry to the free energy of 
the other sphere-like phases of diblock melts 0. This 
scheme has the advantage of providing numerically ex- 
act results using a Fourier basis of a given symmetry 
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FIG. 3: The SCFT phase diagram for mehs composed of (a) 
hnear chains. The dotted hne is where the free energy of the 
BCC phase is equal to that of the A15 phase, while the dashed 
hne is the corresponding boundary for the FCC phase. In the 
phase diagram (b) for the three-generation diblocks shown in 
Fig.n^ the A15 phase is a ground state for a significant range 
of (j> and ■ AH the B blocks are of equal length. Note the 
critical point on the lamellar/disordered boundary. 



and requires significantly less computing power than is 
required to achieve an equivalent precision with a nu- 
merical real-space approach |2^. The efhciency of this 
method allows us to explore the entire phase diagram 
and, because we choose the symmetry of the basis we 
can quickly compare three-dimensional lattices of differ- 
ent symmetries, even in regimes where columnar, lamel- 
lar and gyroid phases are the ground states [i^. Note 



that when considering close-packed lattices, hexagonal 
close-packed and FCC lattices have degenerate free ener- 
gies and so we use the Fm3m basis only. The first result 
of our calculation shows that as ^ the AB interface 
is spherical and as grows, the interface deforms signif- 
icantly into the shape of the Voronoi cell, as shown in 
Fig. 12 This confirms our intuitive discussion about the 
AB interface. 

In order to compute the free energy of all phases within 
0.005% for < 40 we required up to 760 basis func- 
tions. This level of precision allowed us to delineate the 
phase boundaries to within ±0.001 for (p and ±0.01 for 
xN. Our calculation corroborates our Flory-like theory 
and shows that the A15 lattice is the most stable three- 
dimensional lattice for sufhciently large (j), as shown in 
Fig. In addition to the known phases for symmetric 
diblock copolymers, we also computed the boundaries at 
which the free energies of the three candidate sphere-like 
phases crossover. We have reproduced the boundary be- 
tween the BCC and FCC phases reported by Matsen and 
Bates 23|, and have also found the boundaries at which 
the free energy of the A15 lattice becomes lower than that 
of the of FCC and BCC phases, respectively. Thus, we 
see that the A15 lattice is the stable sphere-like phase as 
(j) ^ 1. Note that the boundary between BCC and A15 
terminates a.t (j) = 0.390 and xN — 15.86 because at this 
level of segregation for > 0.390 the Im3m sphere-like 
phase melts. 

Having considered the A15 lattice, it is natural to con- 
sider other periodic arrangements of spheres. We used 
the same SCFT method to compute the free energy of 
sphere-like micelles assembled into the simple cubic (SC) 
and diamond (D) lattice, with Pm3m and Fdim sym- 
metries respectively. For example, at x^ — I'^-O ^-^d 
(f) = 0.26 the free energies per chain in the BCC, FCC, 
A15, SC and D phases are 3.2455, 3.2459, 3.2461, 3.2568 
and 3.2701, respectively in units of ksT. This trend con- 
tinues into the region where sphere-like phases are not 
stable: at = 17-0 and = 0.31 the free energies 
per chain of the A15, BCC, and FCC phases are 3.6107, 
3.6120, and 3.6158, respectively. At this point of the 
phase diagram the free energy per chain of the SC phase 
is 3.6295 and the D phase of spheres becomes unstable. 
The free energies of the BCC, FCC, and A15 phases re- 
main within 1% of each other over the entire region of the 
phase space we explored while the SC and D sphere-like 
phases are never the lowest energy phases. 

Though the A15 lattice will not be observed in neat 
linear diblock melts, our calculation suggests that if the 
sphere-like phases were stabilized it would persist as an 
equilibrium state. One way to stabilize sphere phases 
for more symmetric diblocks is to change the architec- 
ture of the molecule. For instance, if one block of the 
copolymers is branched, the stretching energy combined 
with incompressibility will favor dividing surfaces which 
are curved over a substantial portion of the phase dia- 
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gram. This was predicted by Olmsted and Milner for 
miktoarm block copolymers (l^ and later observed [loj . 
Similarly, calculations by Pickett for multiply branched 
diblocks predict the same generic effect 11]. We have 
implemented a SCFT algorithm to explore the phase di- 
agram of branched diblocks with a linear A-block con- 
nected to n B-blocks, each of which are connected to n 
B-blocks, etc. for g generations so that the final gener- 
ation has {g — l)n B-blocks (the g = 1 system is pure 
A) . We will report elsewhere on the details of our calcu- 
lation (2J] but it closely follows the original analysis of 
Matsen and Schick Q . Because of the time required and 
numerical intensity of the calculation, we have focussed 
first on a 5 = 3 generation diblock with n = 2 blocks at 
each branch as a "proof-of-principle" for both our theo- 
retical argument and our algorithm. As shown in Fig.|3lD, 
we find that sphere-like phases are stable over a greater 
range of (j> and that the A15 is a stable phase at values 
of (j) larger than those for which the BCC phase is sta- 
ble - precisely what we would expect from our geometric 
theory. It is also worth noting that for this region of the 
phase diagram, the gyroid is the equilibrium structure 
for this system, and the double-diamond or hexagonally- 
perforated lamellar phases are only metastable as in lin- 
ear diblocks. 

By employing purely geometric quantities we have ar- 
gued that the A15 lattice should be the lowest energy 
sphere-like micelle phase of diblock copolymers as the 
blocks become more and more symmetric. Through 
a SCFT calculation we have corroborated our theory 
though we found that where the A15 lattice was stable, 
non-spherical columnar or gyroid phases were the ground 
states of the diblock melt. By considering branched di- 
blocks which favor sphere-like phases, we have found via 
SCFT that the A15 is an equilibrium phase of branched 
diblocks with sufficiently large (three) generations. Fur- 
ther work will fully map out the phase diagram of 
branched diblocks as function of molecular geometry and 
topology. 
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